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ABSTRACT 


The purpose of this thesis is to lay the foundation for analyzing structural 
coupling effects for a proposed trolley interface between a ship and a roll-on roll-off 
discharge facility (RRDF). Such a facility could allow heavy cargo transfer at higher sea 
states. Previous studies have analyzed motions assuming that there is no structural 
coupling between the trolley and the RRDF. A mathematical model that incorporates 
structural coupling is developed using the principle of virtual work. In order to assess the 
degree of necessity for the proposed model we conduct a systematic series of numerical 
experiments. In these calculations we model the trolley through a generalized stiffness 
coefficient and assess its influence on RRDF motions. It is shown that modeling of 
structural coupling may be necessary depending on the relative order of magnitude of 


trolley structural rigidity and trolley placement. 
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I. INTRODUCTION 


One of the salient features of the United States Military approach to 21“ Century 
Joint Warfare is Maneuver Warfare. Maneuver warfare seeks to attack adversary 
vulnerabilities from a position of advantage through the synchronized application of 
movement and fires. This is one of the reasons for the creation of the Naval 
Expeditionary Strike Group (NESG). Current NESG comprises of several classes of 
ships such as the LMSR, CAPE-D, CAPE-H and CAPE-T to deliver the necessary 
military equipment and stores to support its mission. In situation where a friendly port is 
available close to the objective, these ships can unload their payload through the port 
facilities. However there are situations where a friendly port is not in close proximity 
resulting in a need to perform at sea transfer of the payload. This involves ship-to-ship 
transfer of payload from these large ships to smaller crafts that subsequently delivers the 


payload to shore. 


Current ship-to-ship transfer involves a ramp interface for heavy equipment such 
as the M1A1 main battle tanks. This technique is limited to operations at sea-state 3 due 
to excessive torsional stresses on the ramp when two main battle tanks are located in the 
proximity of the middle of the ramp. A previous study on a new method of transferring 
cargo based on the use of a trolley-rail interface shows potential advantages in the 
reduction of torsional stresses and weight. This thesis lays the foundation for analyzing 
structural coupling effects for the proposed trolley interface between a ship and a roll-on 
roll-off discharge facility (RRDF), and determines the necessity for modeling of 
structural coupling depending on the relative order of magnitude of trolley structural 


rigidity and trolley placement. 
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Hi. BACKGROUND INFORMATION 


Amphibious operations conducted by the US military involves large amount of 
equipment and stores that needs to be unloaded near to the objective. These equipment 
and stores are usually kept in prepositioned ships at various parts of the world such that 
the response time of the forces can be greatly reduced. However, the current mode of 
operations requires these prepositioned ships to unload their cargo either at nearby port 
facilities, which is the preferred option, or in the event that such ports are not available, 
go right up to within 2 nautical mile of the beach to unload their cargo onto barges or 
small crafts that ferry these equipment and stores ashore. There are two main problems 
with the current mode of operation. First of all, the number of available ports for the 
prepositioned ships to unload their cargo is severely limited. Secondly, the number of 
beaches that would enable these large prepositioned ships to go right near and unload 
their cargo on the beaches is also very limited. The reason for this is that the transfer of 
cargo from the preposition ship to the barge can only be conducted at sea-state three and 
below to prevent excessive torsional stress on the ships’ ramp when unloading heavy 
equipment such as the 68-tonnes M1AI1 main battle tank [Ref 1]. These problems, 
coupled with the fact that such operations results in extended time when the amphibious 
force needs to defend the logistic operations meant that a better way has to be developed. 
This is the main driver for the new Trolley-Rail system proposed by Advanced Design 
Consulting Inc, which aims to enable the transferring of cargo between ships at higher 
sea states. Studies on the response of the proposed system under seaway has yield 
positive results indicating that the system should be able to achieve the required transfer 


rate at higher sea states [Ref 2]. 
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Hil. MODELING 


A. TROLLEY-RAIL MODEL 


The Trolley-Rail system proposed by Advance Design Consulting Inc. essentially 
comprise of two parallel rails that runs from one ship to the other ship. The trolley (or 
platform) traverses along the two rails via attached wheels to transport equipment and 
stores from the prepositioned ships to shuttle ships that ferry these cargo ashore. For the 
model, the other ship is the Roll-on Roll-off Discharge Facility (RRDF). The RRDF is 
made up of a number of interconnected modules which serves as the platform to support 
the rails running from the ship, and to facilitate the movement of the equipment and 
stores to the shuttle ships. A schematic of the Trolley-Rail system is given shown in 
Figure (1) while the three-dimensional model is shown in Figure (2), and the RRDF is 


shown in Figure (3): 





Figure 1. | Schematic of Trolley-Rail System 





-D model of Trolley-Rail System 


3 


Figure 2. 





Schematic of Roll-on Roll-off Discharge Facility (RRDF) 


Figure 3. 


The rails are each modeled as a beam structure. This is a good model as the 
length of each rail is very much greater than its width and thickness. The rails are 
assumed to be pinned at the ship’s end while the other end is free. Each rail is connected 
by a spring, dashpot and actuator to the RRDF to simulate potential active control 
mechanism. The RRDF is connected by a spring and dashpot to the ground to simulate 
the motion of the seaway [Ref 3]. The trolley will effect structural coupling on the two 
rails and this is being modeled as a spring and a dashpot. Two lump masses are also 
attached at the trolley interface with the two rails to account for the cargo loading on the 


trolley. A schematic of the model is shown in Figure 3. 
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> 
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Ship’s End 
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Figure 4. | Model of Trolley-Rail System 


The parameters for the model are given as follows: 

L = length of each Rail 

/ = length between supports along Rail 

L, = length of trolley 

p = distance along the beam from ship end to the center line of the trolley 

Ma = mass contribution of RRDF at the support location of the free end of beam 
k7, = stiffness ratio between Rail # 1 and the RRDF 

C11 = proportional damping ratio between Rail #1 and the RRDF 

kj2 = stiffness ratio between the RRDF and the ground at Rail # | locality 


C12 = proportional damping ratio between the RRDF and the ground at Rail # 1 
locality 


k>, = stiffness ratio between Rail #2 and the RRDF 
C21 = proportional damping ratio between Rail #2 and the RRDF 
kz2 = stiffness ratio between the RRDF and the ground at Rail # 2 locality 


C22 = proportional damping ratio between the RRDF and the ground at Rail # 2 
locality 


v,(x,t) = transverse deflection of Rail #1 

z1(t) = vertical displacement of RRDF at locality of Rail #1 

Z;1(t) = vertical displacement of the sea wave at locality of Rail #1 
v2(x,t) = transverse deflection of Rail #2 

Z2(t) = vertical displacement of RRDF at locality of Rail #2 


Z,2(t) = vertical displacement of the sea wave at locality of Rail #2 


B. DERIVATION OF MATHEMATICAL MODEL 


1. Equation of Motion 


Having developed the schematic of the model, the next step is to derive the 


equations of motion for the model. Lagrange equation was employed for this purpose. 


First of all, the kinetic energy of the system is given by 


T=3I'p 4% a dx +i pa{%) det M, (4 = eee 





Lf a4 Ves ae fees a ss 
Pea ca Piea a 


Working on the premise of small angles, the angular velocity of the RRDF may be 


approximated by, 
ra) = 2, = Z, 
w 
where w is the width of the RRDF. 


Hence, the kinetic energy is given by 





1p yes 2 1{M,)}\,. 2 
beh Af eo 
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The strain energy of the system is given by 





U=— 


é 1 
sf 4) dx+— wk ful t)—z,]) +o hat 


0 





ieeal dx+— aac (I,t)—z,] bokyey 


+k [u(r.d—mv.0} 


For the rails which are modeled as beam structures, the assumed mode method is 


employed to approximate the transverse vibration of each beam, given by 
N 
wnt) = Divide, 
i=l 
N 
V(x 1) = DVB) 
i=l 


where y;(x) is the assumed modes that satisfies the geometric boundary conditions 


of the beam and N is the number of modes that was used in the computation. 


To reduce the mathematical complexity, indicial notations are utilized to represent the 


transverse deflection of the beam. Therefore the transverse deflection of each beam is 


given as 
Vi (4,0) = ya; 
V,(x,t)=W,B, 


where the subscript i represents a summation from | to N. 


Hence, the kinetic energy of the system is represented in the form, 
1 pe ; 
r =a PAly,(a,(O}{y,(a,(O} dx 


“Calpionelly cada 





+ <M, viva} {y, (p)a,(0} |+ <M, wir Boh{y 2,0} | 
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while the potential energy of the system is represented in the form, 
1 L in " 1 2 1 by 
U=5 I, Ely) (xa,()\{y, (xa, (Opde+ ok, [w,Da,(t) - 2] + hat 


+P Ey OBO}; COB O}ar+ rk [VOAO-a] + yhnz 





tk, [vipa O-vi(PBOYy 


Applying the Lagrange equation, 


(ar) ar WU _o 4 
dt \ 0g, 


the generalized coordinates for the model are chosen as, 


~ ‘R 


g=[2, Z, Gh A, ... Ay B, P, sited ae 


The various differentials in the Lagrange equation are given by 


d | OT 1 1 a 
alae (aor) 





a; 


o>) 


(2) {ff ody cow) fer] a, +5 M,C ry, (04, (0 


z ={fi e4[v. cov) te} 5,04 5M vir (8,0 
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or or _ or or 
Oz, Oz, Oa, Of, 





oU 
— =-k,,[y,(~a;,(t)— 2, ]+k,.2, 
Oz, 


= (K,, hs) 2; —kwDa;(t) 


oU 
ae. = (K,, +k,)Z, —kywOBO 
2 


off erly ow; arta 


es ky, [v.(xa, (t) 72 | VY; (/) 


+k, wi(p)a,()-v,(p)B, 0) |v,(P) 


lb elvieweola} ain 


+k, [vi(MB(O-z, |v, 


-k, | wi(p)a,(O-v,(p)B, 0 |v,(P) 


The principle of virtual work is used to determine the generalized force Q;. The 


RRDF and the beams are assumed to have very small movements as follows: 


0Z,, OZ,, Ov,(/,t), and dv,(/,t) at location x =/ 


Ov,(p,t), and dv,(p,t) at location x = p 
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The virtual work done due to dampers resisting force is then given by 
OW" =-c,,2,6z, —¢,,[2,-V,70,0][6z, -dv, LD] 
—CyZ,52, — Cy [Z, —V, (1, |[6z, — dv, 1,0] 


me Cee [v,(p.)-¥,(p,0)][v,(p,0)- dv, (p,0)] 


OW® =[-(e, +¢,,)2, +64, Da, 0] 6z, 
+ | (en + y))Z, + cr DBO | 02, 


+ zy O-cw, Ov, Oa, |6a, 


+ 


l-<, v(Da.0-v,(P)B Oly, (p) |oa, 


+ 


en2wiO-ow, Ov DBO |B, 





+ 


le, {v.(rda,O-v, (PB, O}y,(P) ]5B, 


Hence, the generalize forces due to the dampers are given by 


Oo! = -(c, +0) % +¢,y,0)a,(t) 
oy = (6, HC) 2 t+o,~ (DB (0) 
Of =a yz, +|-awiOv,O+e“,(Py,(p) |@,0-[¢,vi(Dv,(P) |B 


0% =e, (D2, -[¢,vi (Py, (P) |e, 0+] -crvi Ov, O+¢,vi (Pv, (P) |B, 


For the external forces, the virtual work done is given by 


nu, [2acta)(S=82)_ 5 (22-4) 
WwW 





2 2 


Dons 2 
— I, pA 7 Z (t)Ov, (x, t)dx 


ee 
- I, pa rac (t)OV, (x, t)dx 


ae 
~ J, PAT 20M dr |bay, 
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Hence, the generalized forces due to external forces are given by 


2% =-| fasion] in 


24 =-[fieazr cre 20 
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Finally, the equation of motion in generalized coordinates form is given as 


[MN +1C}tg} +1 K]a} = 12) 


whereby the matrices are given by 


M, —-M, O..2N) 
[M] =! —M, M, . 
O22) M wx) 
M B(NXN) 


1 1 
where M, -(4m, +s] 
w 


A A L 1 
M,,=M,, =|, pay, Oy) (x)dx+ 5M vi (PV,(P) 


(A +k.) 0 hiv Oasny Oo) 
0 (ie +k, ) Onxmy ky ¥ i Oaxwny 
[] ~| ky; Oona Oc ay Cjzn) K ben) 
Oca) “hy wv, Oeona Kb en) AuxN) 


A L 

where K,, =| Ely{(xyyi(x)de+ ky, (DvD +k, WP (P) 
A L 

K,, =], Elwin (ade + ky, Ov, O+k, vi Dv,(P) 


Key 7 kW, (Py; (p) 
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(e5 a Gi) 0 CY; Oaxwy 
0 (ey + oi) Dox) 
[C] 5 CW; Day Oona Cas 
Oca CV Oona Cofix) 


where C, =-c,y, Dy) +¢,v,(Dy,(P) 


C,, = envi DvD +¢,vi( DW) 


Cap, =—C VAP, (P) 


My (Zs 259 ) 
M, (24 —2.,) 
{QO} = Pn 
O25 say 
where O, =- i par y,(x)dx 


2. Assumed Modes for Beam 


Oem 
CW Oran 
Cot nny 


An 


ANxN) 


The assumed modes, which satisfy the geometric boundary conditions, need to be 
The 


used to calculate the generalized mass, damping, stiffness and force matrices. 


general solution of the beam is given by [Ref 4] 


y= Acosh £x+ Bsinh Bx+Ccos Bx+ Dsin Bx 
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The boundary conditions are: 




















y=0 
Atx=0 d*y 9 (M = 0) 
dx? 
ay % (M=0) 
Atx=L . 
ae ae = 
dx : ve 
igre) teetac’ 
Py ; pe 7 A=C=0 
[se] -al4 620 A= 
d? 2( Re i 
(2) = p(wsin(ar)-Dsin(p1))=0 
_ sin(BL) D 
~ sinh( BL) 
aq 3 
( 2 J, Sn ee) Dose) =e 
ag cos (BL) D 
cosh (BL) 


sin(ZL) _ cos(BL) 


sinh(BL) cosh( PL) 








= tan( BL) = tanh (BL) 


sin(£,L) 


dy, =C|si 
and y, oi BF) nh) 


sinh ( 3.x) 


In addition, the Eigenfunctions y; utilized in the assumed mode method are 


orthogonal to each other, and therefore the following simplifications can be effected on 


the matrices 


0 (i# /) 
M, =M, = 1 . 
0) pAL+ MVPD) (i=) 
eo ky Oy, O+k,v (Py (Pp) (i# j) 
“| ELAZL+ky(Oy,O+y,(p)y,(P) (i=) 
Bm ky w(Owv (D+k,W, (Py (P) (i# j) 
“| ELIAS L + kw, Ov, D+yv, (Dw (Pp) (i=j) 
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IV. NUMERICAL EXPERIMENTS 


A. METHODOLOGY 


Now that we have developed a mathematical model suitable for structural 


coupling of the trolley/RRDF interface, we need to develop some guidelines with regards 


to its applicability. One main question that we need to answer is whether such a model is 


indeed necessary for this problem. In order to assess this we proceed as follows: 


1. 


First we evaluate the motions of the trolley/RRDF assuming that there is no 
structural coupling. This is achieved in the frequency domain utilizing results 
obtained by WAMIT. WAMIT is a software package that solves the 
hydrodynamic radiation/diffraction problem by a panel discretization of the 
body. It is a 3-dimensional program and as such it can model hydrodynamic 
interactions between multiple bodies. In our case the two bodies are the ship 
and the RRDF. Since we allow for hydrodynamic interactions (in other words 
the motion of each body depends on velocities and accelerations of the 
adjacent body) the equations of motion possess 12 degrees of freedom. Six 
degrees of freedom (surge, sway, yaw, heave, roll, pitch) come from each 


body. 


Based on the computed body motions, we evaluate the relative motions at the 
interface between the two bodies; 1.e., at the trolley end-points. These relative 
motions are subtracted (taking into consideration their relative phase angle) in 


order to obtain the angular displacement of the trolley as well as its twist. 


The next step is to model the structural rigidity of the trolley through a 
generalized stiffness term. The relative motions of the trolley are then used to 
generate the reaction force in the trolley due to its stiffness. This force is then 
added to the hydrodynamic excitation force for the RRDF and the process 
continues. It should be mentioned that such coupling procedure is not 
necessary for the ship since, due to its high displacement, it is unlikely that it 


will be appreciably affected. 
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B. RESULTS 


Results are presented in the next chapter. All results are presented in standard 
polar plots in terms of the RMS value of the indicated response (trolley angle or twist). 
The radial coordinate of the graphs is the significant wave height, while the angular 
coordinate is the wave heading. Fully developed Pierson-Moskowitz and developing or 
decaying Bretchneider wave spectra are used for the calculations. Although the results 
have been obtained for long-crested seas, it is expected that the main conclusions will be 
valid for short-crested seas as well. The trolley stiffness that was introduced in the 
previous section is applied towards the roll motion of the RRDF. This is typical for a 
stern configuration. Similar results will apply for other trolley placements as well. For 
the numerical results that follow, the roll stiffness for the trolley has been normalized 
with respect to the roll stiffness of the RRDF. Naturally, the latter is a function of the 
RRDF displacement and its metacentric height. 
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V. RESULTS 


A. TROLLEY ANGLE AND TWIST SUBJECTED TO _ PIERSON- 
MOSKOWITZ WAVE SPECTRUM 


Figure (5) to (10) shows the RMS value of the trolley angle and trolley twist when 
subjected to the fully developed Pierson-Moskowitz wave spectrum for the cases where 
the normalized trolley stiffness was set to a value of 0, 1 and 2. 180 degrees corresponds 
to head seas while 0 degrees corresponds to following seas. The anti-symmetry in the 
fore and aft responds is due to the placement of the trolley on the aft of the ship. The 
slight anti-symmetry in the Port and Starboard responds is a result of the slight anti- 
symmetry of the shape of the RRDF. It is observed that the trolley angle is not being 
affected by the existence of the trolley stiffness or even if the trolley stiffness is doubled 
in strength. On the other hand, it is observed that the twist experienced by the trolley is 
significantly affected by the introduction of the trolley stiffness coupling term. The effect 


on the trolley twist is amplified further when the normalized trolley stiffness is doubled. 


AMS value of Trolley Angle (deg) 





PM - Rall Stiffness = 0 


Figure 5. RMS Value of Trolley Angle (degrees) subject to Pierce-Moskowitz spectrum and 
zero Roll stiffness 
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RMS value of Trolley Angle {deg} 
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0.2 
270 
PM -Aoll Stiffness = 1 


Figure 6. RMS Value of Trolley Angle (degrees) subject to Pierce-Moskowitz spectrum and 
a Normalized Roll stiffness of 1 
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RAMS value of Trolley Angle {deg} 
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0.2 
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PM - Rall Stiffness = 2 


Figure 7. | RMS Value of Trolley Angle (degrees) subject to Pierce-Moskowitz spectrum and 
a Normalized Roll stiffness of 2 
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RAMS value of Trolley Twist (deg) 





PM - Fall Stiffness = 9 


Figure 8. | RMS Value of Trolley Twist Angle (degrees) subject to Pierce-Moskowitz 
spectrum and zero Roll stiffness 
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RAMS value of Trolley Twist (deg) 





PM - Roll Stiffness = 1 


Figure 9. RMS Value of Trolley Twist Angle (degrees) subject to Pierce-Moskowitz 
spectrum and a Normalized Roll stiffness of 1 
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RAMS value of Trolley Twist (deg) 
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270 


PM - Rall Stiffness = 2 


Figure 10. RMS Value of Trolley Twist Angle (degrees) subject to Pierce-Moskowitz 
spectrum and a Normalized Roll stiffness of 2 
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B. TROLLEY ANGLE AND TWIST SUBJECTED TO BRETSCHNEIDER 
WAVE SPECTRUM 


Figure (11) to (22) shows the RMS value of the trolley angle and trolley twist 
when subjected to the Bretchneider wave spectrum for wave periods of 6 seconds and 12 
seconds. The results are also presented for the three cases where the normalized trolley 
stiffness is set to a value of 0, 1 and 2. In this instance, we observe a similar trend in the 
results with the case under the Pierson-Moskowitz wave spectrum — the trolley angle 
does not seems to be affected by the existence of the trolley stiffness coupling term while 
the trolley twist is substantially affected by the introduction of the trolley stiffness and the 


effects are amplified with an increased in the strength of the coupling term. 


RAMS value of Trolley Angle {deg} 
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0.2 





B6 - Roll Stiffness = 0 


Figure 11. RMS Value of Trolley Angle (degrees) subject to Bretchneider spectrum with a 
wave period of 6 seconds and zero Roll stiffness 
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RAMS value of Trolley Angle {deg} 
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B6 - Roll Stiffness = 1 


Figure 12. RMS Value of Trolley Angle (degrees) subject to Bretchneider spectrum with a 
wave period of 6 seconds and a Normalized Roll stiffness of 1 
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RAMS value of Trolley Angle {deg} 
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B6 - Rall Stiffness = 2 


Figure 13. RMS Value of Trolley Angle (degrees) subject to Bretchneider spectrum with a 
wave period of 6 seconds and a Normalized Roll stiffness of 2 
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RAMS value of Trolley Twist (deg) 





B6 - Roll Stiffness = 0 


Figure 14. RMS Value of Trolley Twist Angle (degrees) subject to Bretchneider spectrum 
with a wave period of 6 seconds and zero Roll stiffness 


RAMS value of Trolley Twist (deg) 











Bé - Roll Stiffness = 4 


Figure 15. RMS Value of Trolley Twist Angle (degrees) subject to Bretchneider spectrum 
with a wave period of 6 seconds and a Normalized Roll stiffness of 1 
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RAMS value of Trolley Twist (deg) 
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B6 - Rall Stiffness = 2 


Figure 16. RMS Value of Trolley Twist Angle (degrees) subject to Bretchneider spectrum 
with a wave period of 6 seconds and a Normalized Roll stiffness of 2 
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RAMS value of Trolley Angle {deg} 
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B12 - Rall Stiffness = 0 


Figure 17. RMS Value of Trolley Angle (degrees) subject to Bretchneider spectrum with a 
wave period of 12 seconds and zero Roll stiffness 
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RAMS value of Trolley Angle {deg} 
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Figure 18. RMS Value of Trolley Angle (degrees) subject to Bretchneider spectrum with a 
wave period of 12 seconds and a Normalized Roll stiffness of 1 
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RAMS value of Trolley Angle {deg} 
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Figure 19. RMS Value of Trolley Angle (degrees) subject to Bretchneider spectrum with a 
wave period of 12 seconds and a Normalized Roll stiffness of 2 
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RAMS value of Trolley Twist (deg) 





Bi2 - Rall Stiffness = 0 


Figure 20. RMS Value of Trolley Twist Angle (degrees) subject to Bretchneider spectrum 
with a wave period of 12 seconds and zero Roll stiffness 


RAMS value of Trolley Twist (deg) 





B12 - Roll Stiffness = 1 


Figure 21. RMS Value of Trolley Twist Angle (degrees) subject to Bretchneider spectrum 
with a wave period of 12 seconds and a Normalized Roll stiffness of 1 
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RAMS value of Trolley Twist (deg) 
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Figure 22. RMS Value of Trolley Twist Angle (degrees) subject to Bretchneider spectrum 
with a wave period of 12 seconds and a Normalized Roll stiffness of 2 
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C. INFLUENCE OF TROLLEY STIFFNESS ON TROLLEY TWIST 


The observations from the previous sections are better illustrated through Figure 
(23) to (28), which shows the comparison of the responds of the trolley angle and trolley 
twist, with and without the trolley stiffness coupling term. The almost symmetrical 
response of Port and Starboard direction is expected with the slight variation being 
accounted for by the slight anti-symmetry of the shape of the RRDF. It can clearly be 
seen that there is substantial differences in the responds of the trolley twist when the 
stiffness coupling term is introduced. This responds is magnified when the coupling term 
is doubled. In addition, it is observed that the quartering seas results in higher trolley 


twist. 


Diff RMS Trolley Twist {deg} 





270 


PM - Rall Stiffness - 0 vs 1 


Figure 23. Difference in average vertical trolley twist angle for a Normalized Roll stiffness 
of 0 and 1 (subjected to Pierce-Moskowitz wave spectrum) 
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Diff RMS Trolley Twist {deg} 





PM - Rall Stiffness - 0 vs 2 


Figure 24. Difference in average vertical trolley twist angle for a Normalized Roll stiffness 
of 0 and 2 (subjected to Pierce-Moskowitz wave spectrum) 
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Figure 25. Difference in average vertical trolley twist angle for a Normalized Roll stiffness 
of 0 and 1 (subjected to Bretchneider spectrum with a wave period of 6 seconds) 
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Diff RMS Trolley Twist {deg} 
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B6 - Rall Stiffness - 0 vs 2 


Figure 26. Difference in average vertical trolley twist angle for a Normalized Roll stiffness 
of 0 and 2 (subjected to Bretchneider spectrum with a wave period of 6 seconds) 
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Diff RMS Trolley Twist {deg} 
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Figure 27. Difference in average vertical trolley twist angle for a Normalized Roll stiffness 
of 0 and 1 (subjected to Bretchneider spectrum with a wave period of 12 seconds) 
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Diff RMS Trolley Twist {deg} 
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Bi2 - Roll Stiffness - 0 vs 2 


Figure 28. Difference in average vertical trolley twist angle for a Normalized Roll stiffness 
of 0 and 2 (subjected to Bretchneider spectrum with a wave period of 12 seconds) 
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VI. CONCLUSIONS AND RECOMMENDATIONS 


A. CONCLUSIONS 


This thesis has laid the foundation for the analyzing of structural coupling effects 
for the proposed trolley interface between a ship and a Roll-on Roll-off Discharge 
Facility (RRDF) through the development of a mathematical model utilizing Assumed 
Modes and Virtual Work techniques, which incorporates the structural coupling effects of 
the trolley to the RRDF. Through a systematic series of numerical experiments, it was 
also determined that the trolley angle is not affected by the structural coupling effects, in 


contrast with the trolley twist which is substantially affected by structural coupling. 


B. RECOMMENDATIONS FOR FUTURE WORK 


1. Side Trolley Placement 


We have conducted the analysis for stern placement of the new trolley-rail 
system. However there are cases where the trolley-rail system may be placed on the side 
of the ship for unloading operations. In this scenario, the hydrodynamic coupling effects 
of the ship and the RRDF will be significantly different. A similar analysis should be 
conducted for the side trolley placement to determine the relative impact of the structural 


coupling effects. 


2. Structural Coupling for Trolley Twist 


Given the knowledge that the trolley twist is substantially affected by structural 
coupling between the trolley and RRDF, we can now proceed to analyze the structural 
coupling effects utilizing the mathematical model developed, so that appropriate 
measures can be taken to ensure that cargo transfer can be conducted at sea up to Sea 


state 5. 
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APPENDIX A 


MATLAB CODE 


% Motions calculation for CapeD/RRDF’/trolley interface. 
% Long-crested Bretchneider seas. 
% 
i_seaway=input('Enter 0 for PM or | for Bretchneider spectrum = '); 
if i_ seaway == 1 

T_m =input('Modal Period (sec) = '); 

omega _ m=2*pi/T_m; 
end 
i_stiffness=input('Enter 3 for heave stiffness, 4 for roll, or 5 for pitch = '); 
k_stiffness=input('Enter appropriate normalized stiffness = '); 
% 
x_arm_CAPED_1=-330; 
x_arm_CAPED 2=-330; 
y_arm_CAPED _1=-10; 
y_arm_CAPED 2=+10; 
x_arm_RRDF_ 1=20; 
x_arm_ RRDF 2=20; 
y_arm_ RRDF _ 1=-10; 
y_arm RRDF 2=+10; 
trolley_length=100; 
% 
% Load processed WAMIT data file 
% 
load CAPED 
NH=25; 
i NH=15; 
beta_incr=i_NH; 
imag_unit=sqrt(-1); 
% 
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% Set the frequencies vector (0.3 to 2.5 rad/sec) 
% 
NF=(size(frequencies)); 
NF=NF(1,1); 
w=frequencies; 
% 
% Set the headings vector (0 to 360 degrees in i_ NH deg. increments) 
% 
for i=1:NH, 
heading(i)=i_ NH*(i-1); 
end 
% 
% Set added mass and damping matrices and forcing vector 
% 
for i=1:NF, 
i_string=num?str(i); 
A(.,:,1)=eval(streat(‘admassfreq',i_string)); 
BC:,:,1)=eval(strceat(‘addumpingfreq',1_string)); 
for j=1:NH, 


j_string=num2str(j); 


j_string=num2str(1); 
end 
F(:,1,j)=eval(strcat('forcfreq',i_string,'head',j_string)); 
end 
end 
% 
% Random wave calculations 
% PM spectrum - weighted by the spreading function 
% 
% Loop on significant wave height 
iHS=0; 


for HS=0.5:0.5:30, 
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HS; 
iHS=iHS+1; 
if i_ seaway ==0 
omega_m=0.4*sqrt(32.2/HS); 
end 
A_s=(1.25/4)*(omega_m*4)*(HS%2); 
B_s=1.25*omega_ m4; 
S_main=(A_s./w.4%5).*exp(-B_s./w.%4); 
% 
% Loop on sea direction 
% 
ibeta=0; 
for beta=0:beta_incr:360, 
beta_set=beta; 
ibeta=ibetat1; 
S=S_main; 
heading _single=heading(ibeta); 
% 


% Frequency domain response is given by x= [-w42*(A+mass)+i*w*B+C]\(-1)*F 





% Zeroth order calculation (no trolley/RRDF structural coupling) 
% 
for i=1:NF, 
w_single=w(i); 
A_bar=-w_single*w_single*(A(:,:,i))+imag_unit*w_single*B(:,:,1)+C; 
F_bar=F(:,i,ibeta); 
x=inv(A_bar)*F_bar; 
% 
% Ship motions 
% 
surge0_ CAPED(i) = x(1); 
sway0_ CAPED(i) = x(2); 
heaveO_ CAPED(i) = x(3); 


roll0. CAPED(i) = x(4); 
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pitchO. CAPED() = x(5); 
yaw0_CAPED(i) =x(6); 

% 

% RRDF motions 

% 

surge0_ RRDF(i) = x(7); 

swayO RRDF(1) = x(8); 

heaveO_ RRDF(i) = x(9); 

roll0 RRDF(i) = x(10); 

pitchO RRDF(i) = x(11); 

yaw0 RRDF(i) =x(12); 

% 

% Motions/ at trolley end - CAPED side 
% 

m0_t_ CAPED v_1(i)=heave0Q_CAPED(1)- 


x_arm_CAPED _1*pitchO CAPED(i)... 


-y_arm_CAPED_1*roll0_ CAPED(); 
m0_t_ CAPED v_2(i)=heave0_CAPED(1)- 


x_arm_CAPED 2*pitchO CAPED(i)... 


m0 _t CAPED h_ 1(i)=abs(swayO CAPED(i))+x_arm_ CAPED_1*abs(yaw0_CAPED(i)) 


> 


m0 _t CAPED h 2(i)=abs(swayO CAPED(i))+x_arm CAPED 2*abs(yaw0_ CAPED(i)) 


> 


m0 t RRDF h_1(i)=abs(swayO_RRDF(i))+x_arm_ RRDF_1*abs(yaw0_RRDF(i)); 


-y_arm_CAPED 2*roll0_ CAPED(i); 


% 
% Motions at trolley end - RRDF side 
% 


m0_t RRDF_v_1(i)=heave0_RRDF(i)-x_arm_ RRDF_1*pitchO RRDF(i)... 


-y_arm RRDF_ 1*roll0 RRDF(); 


m0_t RRDF_v_2(i)=heave0_RRDF(i)-x_arm_ RRDF_ 2*pitchO RRDF(i)... 


-y_ arm RRDF 2*roll0 RRDF(i); 
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m0 _t RRDF h 2(i)=abs(swayO_RRDF(i))+x_arm_ RRDF 2*abs(yaw0_RRDF(i)); 
% 
% Trolley relative vertical angular displacement 
% 


m0 trolley angle v_1(i)=(abs(heave0_ CAPED(i))- 
x_arm_CAPED _1*abs(pitchO_CAPED(i))... 


-y_arm CAPED _1*abs(roll0_ CAPED(i))-(abs(heave0O_ RRDF(i))- 
x_arm_RRDF _1*abs(pitchO_RRDF(i))... 


-y_arm_RRDF _1*abs(roll0_ RRDF(i))))/trolley_length; 


m0 trolley angle v_2(i)=(abs(heave0_ CAPED(i))- 
x_arm_CAPED 2*abs(pitchO_CAPED(i))... 


-y_arm CAPED 2*abs(roll0_ CAPED(1))-(abs(heave0O_ RRDF(i))- 
x_arm_RRDF 2*abs(pitchO RRDF())... 


-y_arm_RRDF 2*abs(roll0_ RRDF(i))))/trolley_length; 
% 
% Trolley relative horizontal displacement 
% 
m0_trolley_distance_h_1(i)=(m0_t_CAPED h 1(i)-m0_t RRDF_h_1(i)); 
m0_trolley_distance_h 2(i)=(m0_t_ CAPED h 2(i)-m0_ t RRDF_h_ 2(i)); 
% 
% Calculate transmission force due to trolley end-points motion differential 
% 
if i_stiffness==3 
f diff(ij=C(9,9)*k_stiffness*(m0_t RRDF_v_1(i)-m0_t RRDF_v_2(i)); 
heave_trolley(i)=f_diff(); 
roll_trolley(i)=0; 








pitch_trolley(i)=0; 
end 
if i_stiffness==4 


f diff(~ij=C(10,10)*k_stiffness*abs((m0_trolley_angle v_1(i)- 
m0 trolley angle v_2(1))); 








heave_trolley(i)=0; 
roll _trolley()=f_diff(i); 
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pitch_trolley(i)=0; 
end 
if i_stiffness== 


f diff(ij=C(11,11)*k_stiffness*abs((m0_trolley_angle v_1(i)- 
m0 trolley angle v_2(1))); 








heave_trolley(i)=0; 
roll_trolley(i)=0; 
pitch_trolley(i)=f_diff(i); 
end 
end 
% 
% First-order correction 
% Calculate motions again incorporating trolley stiffness correction 
% 
for i=1:NF, 
w_single=w(i); 
A_bar=-w_single*w_single*(A(:,:,1))t+imag_unit*w_single*B(:,:,1)+C; 
f trolley=f_diff(i); 
m_trolley=f_diff(i)*20; 


F_bar=F(:,i,ibeta)+[0;0;0;0;0;0;0;0;heave_trolley(i);- 
roll_trolley(i);pitch_trolley(i);0]; 


x=Inv(A_bar)*F_bar; 

% 

% Ship motions 

% 

surge CAPED(i) = x(1); 
sway_CAPED(i) = x(2); 
heave CAPED(i) = x(3); 
roll CAPED(i) = x(4); 
pitch CAPED(i) = x(5); 
yaw_CAPED(i) = x(6); 
% 

% RRDF motions 
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% 

surge RRDF(i) = x(7); 

sway_RRDF(i) = x(8); 

heave RRDF(i) = x(9); 

roll RRDF(i) = x(10); 

pitch RRDF(i) = x(11); 

yaw RRDF(i) =x(12); 

% 

% Motions/ at trolley end - CAPED side 

% 

m_t CAPED v_ 1(i)=heave CAPED(i)-x_arm CAPED_1*pitch CAPED(i)... 
-y_arm_CAPED_1*roll_ CAPED(i); 

m_t CAPED v 2(i)=heave CAPED(i)-x_arm CAPED 2*pitch CAPED(i)... 
-y_arm_CAPED 2*roll_ CAPED(i); 


m_t CAPED h 1(i)=abs(sway_ CAPED(i))+x_arm CAPED_1*abs(yaw_CAPED(i)); 





m_t CAPED h 2(i)=abs(sway_CAPED(i))+x_arm CAPED 2*abs(yaw_CAPED(i)); 
% 
% Motions at trolley end - RRDF side 
% 
m_t RRDF v_1(i)=heave RRDF(i)-x_arm_ RRDF_1*pitch RRDF(i)... 
-y_arm RRDF_ 1*roll RRDF(i); 
m_t RRDF v_2(i)=heave RRDF(i)-x_arm_ RRDF_2*pitch RRDF(i)... 
-y_arm RRDF 2*roll RRDF(i); 
m_t RRDF h_1(i)=abs(sway_RRDF(i))+x_arm_RRDF_1*abs(yaw_RRDF(i)); 
m_t RRDF h 2(i)=abs(sway_RRDF(i))+x_arm_RRDF_2*abs(yaw_RRDF(i)); 
% 
% Trolley relative vertical angular displacement 
% 
m_trolley angle v_1(i)=(abs(heave CAPED(i))- 
x_arm_CAPED_1*abs(pitch_CAPED(i))... 
-y_arm CAPED _1*abs(roll_ CAPED(i))-(abs(heave_RRDF(i))- 
x_arm_RRDF_1*abs(pitch RRDF(i))... 














a3 


-y_arm_RRDF _1*abs(roll_ RRDF(1))))/trolley_length; 


m_trolley angle v_2(i)=(abs(heave CAPED(i))- 
x_arm_CAPED 2*abs(pitch_CAPED(i))... 


-y_arm_ CAPED 2*abs(roll_ CAPED(i))-(abs(heave_RRDF(i))- 
x_arm_RRDF _ 2*abs(pitch RRDF(i))... 


-y_arm_RRDF 2*abs(roll_ RRDF(1))))/trolley_length; 


m_trolley angle v_twist(i)=m_trolley angle v_1(i)-m trolley angle v_2(i); 











m_trolley angle v_average(i)=0.5*(m_trolley_angle v_1(i)+m trolley angle v_2(a)); 

% 
% Trolley relative horizontal displacement 
% 
m_trolley distance h 1(i)=(m_t_ CAPED h 1(i)-m t RRDF h 1(i)); 
m_trolley distance h 2(i)=(m_t CAPED h 2(i)-m t RRDF h 2(i)); 

end 

% 

% Define response spectra 

% 

S_ surge CAPED = ((abs(surge CAPED)).%2).*S'; 

S_ surge RRDF = ((abs(surge RRDF)).%2).*S'; 

S_ heave CAPED = ((abs(heave CAPED)).‘2).*S'; 

S_ heave RRDF = ((abs(heave RRDF)).‘2).*S'; 

S_sway_CAPED = ((abs(sway_CAPED)).‘2).*S'; 

S sway RRDF = ((abs(sway_RRDF)).%2).*S'; 

S_roll CAPED = ((abs(roll CAPED)).‘2).*S'; 

S roll RRDF = ((abs(roll_ RRDF)).%2).*S'; 

S_pitch CAPED = ((abs(pitch CAPED)).%2).*S'; 

S_ pitch RRDF = ((abs(pitch_RRDF)).%2).*S'; 

S_yaw_CAPED = ((abs(yaw_CAPED)).‘2).*S'; 

S yaw _RRDF = ((abs(yaw_RRDF)).%2).*S'; 

S_m_t CAPED v_1=((abs(m_ t CAPED v_1)).%2).*S 

S_m_t CAPED v 2=((abs(m_ t CAPED v_2)).%2).*S'; 

S_m t CAPED h 1=((abs(m_t CAPED h_ 1)).‘%2).*S'; 

















54 


S_m_t CAPED h 2=((abs(m_ t CAPED h 2)).%2).*S'; 
S_m_t RRDF v_1 =((abs(m_t RRDF v_1)).%2).*S5 

S_m_t RRDF v 2 =((abs(m_t RRDF v _ 2)).%2).*S5 

S_m_t RRDF h 1 =((abs(m_t RRDF h_1)).%2).*S5 

t RRDF h 2 =((abs(m_t RRDF h 2)).42).*S'; 

trolley angle v_1 =((abs(m_trolley_angle v_1)).42).*S'; 




















trolley angle v_2 =((abs(m_trolley_angle v_2)).%2).*S'; 








m 
m 
m 
m 


trolley angle v_twist =((abs(m_trolley_angle v_twist)).42).*S'; 








S 
S 
S 
Ss 
S_m trolley angle v_ average = ((abs(m_trolley_angle v_average)).2).*S'; 
S_m trolley distance _h_1 = ((abs(m_trolley_distance_h_1)).%2).*S'; 

S_m trolley distance h 2 =((abs(m trolley distance h 2)).%2).*S'; 

% 

% Integral initializations 

% 

S_ surge CAPED i=0; 

S_ surge RRDF i =0; 

S_heave CAPED i=0; 

S heave RRDF i =0; 

S_sway_CAPED i =0; 

S_sway_RRDF i =0; 

S_roll CAPED i =0; 

S_roll RRDF i =0; 

S_pitch_ CAPED i=0; 

S_pitch RRDF i =0; 

S_yaw_CAPED i =0; 

S yaw RRDF i =0; 

S_m t CAPED v_1 i=0; 

S_m t CAPED h 1 i=0; 

S mt RRDF v 1 i =0; 

S mt RRDF h 1 i =0; 

S_m _trolley_angle v_1 i=0; 





























S_m_trolley_distance_h_1 i=0; 


S_m_t CAPED v 2 i=0; 
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S_m_t CAPED h 2 i1=0; 
Smt RRDF v 2 i =0; 
S mt RRDF h 2 i =0; 
S_m _trolley_angle v 2 i=0; 











S_m_trolley_distance_h 2 i=0; 


S_m trolley angle v_twist_i=0; 





S_m_trolley_angle v_average i= 0; 

% 

% Integral S(w)*|RAO(w)|*2 

% 

for i=2:1:NF, 
avg value=0.5*(S_surge CAPED(i)+S_ surge CAPED(i-1)); 
S_ surge CAPED i=S_ surge CAPED itavg value*abs(w(i-1)-w()); 
avg value=0.5*(S_surge RRDF(i)+S_ surge RRDF(i-1)); 
S_ surge RRDF i =S_ surge RRDF itavg value*abs(w(i-1)-w(i)); 
avg value=0.5*(S_ heave CAPED(i)+S_ heave CAPED(i-1)); 
S heave CAPED i=S heave CAPED itavg value*abs(w(i-1)-w(i)); 
avg value=0.5*(S_ heave RRDF(i)+S_ heave RRDF(i-1)); 
S heave RRDF i =S_ heave RRDF itavg value*abs(w(i-1)-w(i)); 
avg value=0.5*(S_sway CAPED(i)+S_sway_CAPED(i-1)); 
S_ sway CAPED i =S_ sway _CAPED itavg value*abs(w(i-1)-w(i)); 
avg value=0.5*(S_ sway RRDF(i)+S_sway_RRDF(i-1)); 
S_ sway RRDF i =S_ sway _RRDF it+avg value*abs(w(i-1)-w(i)); 
avg value=0.5*(S_roll CAPED(i)+S_roll CAPED(-1)); 
S_ roll CAPED i =S roll CAPED itavg value*abs(w(i-1)-w(i)); 
avg value=0.5*(S_roll RRDF(i)+S_roll RRDF(i-1)); 
S roll RRDF 1 =S roll RRDF itavg_ value*abs(w(i-1)-w(i)); 
avg value=0.5*(S_ pitch CAPED(i)+S_pitch CAPED(i-1)); 
S pitch CAPED i=S pitch CAPED itavg value*abs(w(i-1)-w(i)); 
avg value=0.5*(S_pitch RRDF(1)+S_pitch RRDF(i-1)); 
S pitch RRDF i =S_ pitch RRDF i+avg value*abs(w(i-1)-w(i)); 
avg value=0.5*(S_yaw_CAPED(i)+S_ yaw CAPED(i-1)); 


S_yaw_CAPED i =S_yaw_CAPED itavg_value*abs(w(i-1)-w(1)); 
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avg value=0.5*(S_yaw_RRDF(i)+S_ yaw _RRDF(i-1)); 

S_yaw_RRDF i =S_ yaw RRDF itavg value*abs(w(i-1)-w(i)); 

% 

% Calculations for point "1" 

% 

avg value=0.5*(S_m t CAPED v_ 1(1)+S m t CAPED v_ l(i-1)); 

S mt CAPED v_11=S mt _CAPED v _1 itavg value*abs(w(-1)-w(i)); 
avg value=0.5*(S_m t CAPED h 1(1)+S m t CAPED h 1(i-1)); 

S mt CAPED h 1 1=S mt CAPED h 1 itavg value*abs(w(-1)-w(i)); 
avg _value=0.5*(S_m t RRDF v_1()+S m t RRDF v_l(i-1)); 

S mt RRDF v_1i=S mt RRDF v_1 ttavg value*abs(w(i-1)-w()); 
avg value=0.5*(S_m t RRDF h 1(i)+S mt RRDF h 1(i-1)); 

Smt RRDF h 1 i=S mt RRDF h 1 ttavg value*abs(w(i-1)-w()); 









































avg value=0.5*(S_m trolley angle v_1(i)+S m trolley angle v_1(i-1)); 





S_m_ trolley angle v_1i=S_m trolley angle v_1 1+... 





avg value*abs(w(i-1)-w(i)); 


avg value=0.5*(S_m trolley distance h 1(i)+S_ m trolley distance h_1(1-1)); 





S_m trolley distance _h_1i=S_m trolley distance h_1 i+... 





avg value*abs(w(i-1)-w(i)); 
% 
% Calculations for point "2" 
% 
avg value=0.5*(S_m t CAPED v 2(1)+S m t CAPED v 2(i-1)); 
S mt CAPED v 2 1=S mt CAPED v 2 itavg value*abs(w(1-1)-w(i)); 
avg value=0.5*(S_m t CAPED h 2(1)+S m t CAPED h 2(i-1)); 
S mt CAPED h 2 1=S mt CAPED h 2 itavg value*abs(w(1-1)-w(i)); 
avg _value=0.5*(S_m t RRDF v_ 2(i)+S m t RRDF v_2(i-1)); 
S mt RRDF v 2 i1=S mt RRDF v 2 ttavg value*abs(w(i-1)-w()); 
avg value=0.5*(S_m t RRDF h 2(i)+S mt RRDF h 2(i-1)); 
S mt RRDF h 2 i=S mt RRDF h 2 ttavg value*abs(w(1-1)-w()); 






































avg value=0.5*(S_m trolley angle v_2(i)+S m trolley angle v_2(i-1)); 





S_m_trolley_angle v 2 i=S_m trolley angle v 2 i+... 





avg value*abs(w(i-1)-w(i)); 
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avg value=0.5*(S_m trolley distance h 2(i)+S_ m trolley distance h 2(1-1)); 





S_m trolley distance h 2 i=S_m trolley distance _h 2 i+... 





avg value*abs(w(i-1)-w(i)); 
% 
avg value=0.5*(S_m trolley angle v_twist(i)+S m_ trolley angle v_twist(i-1)); 








S_m _trolley_angle v_twist_i=S_m trolley angle v_twist_it... 





avg value*abs(w(i-1)-w(i)); 


avg value=0.5*(S_m trolley angle v_average(i)+S m_ trolley angle v_average(i-1)); 








S_m trolley angle v average 1=S _m trolley angle v_ average it... 








avg value*abs(w(i-1)-w(i)); 
end 
% 
% RMS values 
% 
RMS surge CAPED(ibeta,iHS)= sqrt(S_ surge CAPED_ 1); 
RMS surge RRDF(ibeta,iHS) = sqrt(S_surge RRDF _1); 
RMS _heave_CAPED(ibeta,iHS)= sqrt(S heave CAPED_ 1); 
RMS heave RRDF(ibeta,iHS) = sqrt(S heave RRDF 1); 
RMS _sway_CAPED(ibeta,iHS) = sqrt(S_ sway _CAPED 1); 
RMS _ sway RRDF(ibeta,iHS) = sqrt(S sway RRDF _ i); 
RMS roll CAPED(ibeta,iHS) = sqrt(S_roll_ CAPED_1); 
RMS roll RRDF(ibeta,iHS) = sqrt(S_roll_ RRDF_1i); 
RMS pitch CAPED(ibeta,iHS)= sqrt(S_pitch CAPED_1); 
RMS pitch RRDF(ibeta,iHS) = sqrt(S_ pitch RRDF_1); 
RMS_yaw_CAPED(ibeta,iHS) = sqrt(S_yaw_CAPED_ 1); 
RMS_yaw_RRDF(ibeta,iHS) =sqrt(S yaw _RRDF _ 1); 
RMS _m_t_CAPED v_l(ibeta,iHS) =sqrt(S_m_t CAPED v_1 1); 
RMS _m_t_CAPED h l(ibeta,iHS) =sqrt(S_m_t_CAPED h 1 1); 
RMS m t RRDF v_ l(ibeta,iHS) =sqrt(S m t RRDF v_1 1); 
RMS m t RRDF h l(ibeta,iHS) =sqrt(S m t RRDF h 1 1); 
RMS m trolley_angle v_1(ibeta,iHS) = sqrt(S_m_trolley_angle v_1 1); 
RMS m trolley_distance_h_1(ibeta,iHS) = sqrt(S_m_trolley_distance_h_1_1); 
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RMS _m_t_CAPED v_ 2(ibeta,iHS) = sqrt(S_m_t CAPED v 2 1); 
RMS _m_t_CAPED h 2(ibeta,iHS) = sqrt(S_m_t CAPED h 2 1); 
RMS m t RRDF v 2(ibeta,iHS) =sqrt(S m t RRDF v 2 1); 
RMS m t RRDF h 2(ibeta,iHS) =sqrt(S m t RRDF h 2 1); 














RMS m trolley_angle v_2(ibeta,iHS) = sqrt(S_m_trolley_angle v 2 1); 
RMS m trolley_distance_h_ 2(ibeta,iHS) = sqrt(S_m_trolley_distance_h_2_1); 








RMS m trolley_angle v_twist(ibeta,iHS) = sqrt(S_m_trolley_angle_v_twist_1); 








RMS m trolley_angle v_average(ibeta,iHS) 
sqrt(S_m_trolley_angle_v_average_1); 








end 
end 
% 
figure(1) 
[th,r]=meshgrid((0:beta_incr:360)*pi/180,0.5:0.5:30); 
[X, Y |=pol2cart(th,r); 
h=polar(th,r);delete(h); 
hold on 
c_p=[0.0:0.1:2.0]; 
contour(X',Y',\RMS_m trolley angle v_average*180/pi,c_p),caxis([0 2]),colorbar 
title(RMS value of Trolley Angle (deg)') 
% 
figure(2) 
[th,r]=meshgrid((0:beta_incr:360)*pi/180,0.5:0.5:30); 
[X, Y |=pol2cart(th,r); 
h=polar(th,r);delete(h); 
hold on 
c_p=[0.0:0.1:2.0]; 
contour(X',Y'.\RMS_m trolley angle v_twist*180/pi,c_p),caxis([0 2]),colorbar 
title((RMS value of Trolley Twist (deg)') 
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